Bond distortion effects and electric orders in spiral multiferroic magnets 
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We study in this paper bond distortion effect on electric polarization in spiral multiferroic magnets 
based on cluster and chain models. The bond distortion break inversion symmetry and modify the 
d-p hybridization. Consequently, it will affect electric polarization which can be divided into spin- 
current part and lattice-mediated part. The spin-current polarization can be written in terms of 
Si,j x (S x e,) and the lattice-mediated polarization exists only when the M-O-M bond is distorted. 
The electric polarization for three-atom M-O-M and four-atom M-O2-M clusters is calculated. We 
also study possible electric ordering in three kinds of chains made of different clusters. We apply 
our theory to multiferroics cuprates and find that the results are in agreement with experimental 
observations. 



I. INTRODUCTION 

Magnetoelectric multiferroics, where magnetic and 
electric ordering coexisting in a single compound, has at- 
tracted much interest in the last decades due to their po- 
tential for novel physics and technological applications 1 . 
However, the initial observations have shown that ma- 
terials with coexisting ferroelectric and magnetic orders 
are rare in nature and the observed magnetoelectric cou- 
pling was rather weak 2 . Therefore, the search for new 
gigantic magnetoelectric multiferroic materials is a great 
challenge. As an important progress in multiferroics, 
Kimura et ai. 3 found giant magnetoelectric effects in the 
perovskite manganite TbMnC>3 in 2003. The ferroelec- 
tric transition at T c = 28 K in TbMnC>3 was ascribed to 
the emergence of a nontrivial magnetic order, i.e., with 
a spiral spin structure^. This finding revived the inter- 
ests in multiferroic behavior and lead to the discovery 
of a relatively large number of spin-driven multiferroic 
materials^—. 

One of the fundamental issues in multiferroic physics 
is to clarify the origin of the magnetoelectric coupling in 
such new materials. Currently, although several mod- 
els have been proposed both microscopically^— and 
phenomenologicallyi&, the microscopic mechanism of the 
magnetoelectric coupling is still a controversial and un- 
resolved topic. One of the most prevailing mechanisms 
is the spin-current model proposed by Katsura et al. s 
for helical magnets. They studied a three-atom M-O-M 
(M denotes a transition metal ion and O is for an oxy- 
gen ion) cluster where O locates in the inversion center of 
the cluster, say, the mid point between two M ions. They 
pointed out that spontaneous electric polarization in the 
spiral magnetic order arises from a distortion of electronic 
density even without ionic or atomic displacements. And 
the polarization can be expressed in the form 



P oc ei,j x (ei x &j), 



(1) 



where e^j denotes the bond direction connecting the two 
neighbor spin moments Si and Sj along directions e"i 
and ej on the sites i and j respectively. The electric 



polarization is closely associated with the nonzero spin 
current Si x Sj between two neighboring noncollinear 
spins. Furthermore, Jia et al£^ presented a more so- 
phisticated calculation on the spin-current model with 
realistic considerations. Their results showed that the 
spin-current model is able to explain many experimen- 
tal data at least semiquantitatively. On the other hand, 
Sergienko and Dagotto^ argued that the ferroelectric po- 
larization is closely related to the inverse Dzyaloshinskii- 
Moriya(DM) coupling which breaks inversion symme- 
try and obtained similar expression for electric polariza- 
tion as in ([1]). Interestingly, the spin-orbit interaction 
plays a crucial role in both theories. The role of spin- 
orbit interaction was also confirmed by Mostovoy based 
on the phenomenological theory 18 . Experimentally, the 
spin-current mechanism was demonstrated directly by 
controlling the spin-helicity vector in TbMn03 with an 
external electric field 19 . The spin-current mechanism 
has also been applied to several other spiral multiferroic 
materials, such as DyMnO^, MnWO^, NiaV^Og^, 
CoCr 2 OA LiCu 2 O^^ am i LiCuVO^- 2 ^. Beside 
these successful applications, it should be noted that the 
magnitude of electric polarization is very sensitive to the 
location of Fermi energy in existing spin-current models. 
When the Fermi energy locates between two hybridized 
levels, which refers to "single hole" in Ref* 8 -, P ex V/A, 
where V is the hybridization between transition metal d- 
orbitals and oxygen p-orbitals and A is the energy differ- 
ence between d-orbitals and p-orbitals. Whereas, when 
the Fermi energy locates outside two hybridized levels, 
which refers to "double hole" in Ref. 8 , P oc (V/A) 3 , 
which is significantly different from "single hole" situ- 
ation. 

In real multiferroics, their fascinating properties are 
attributed to the competition among charge, spin, or- 
bital, and lattice degrees of freedom. Actually, the 
bond-distortion is inevitably present in realistic transi- 
tion metal oxides due to structure distortio n 30 ' 31 . Thus, 
the inversion symmetry breaks and the oxygen atom is 
away from the mid point between two M atoms. There- 
fore, the bond-distortion effects should be taken into ac- 
count in the calculation for the electric polarization. In 



existing spin-current model^— , the bond-distortion ef- 
fect is overlooked. As a result that we shall show in this 
paper, the calculated electric polarization is underesti- 
mated. A recent theoretical study on a toy model for 
the ferroelectricity of a two-dimensional cluster— sug- 
gests that the bond-bending may be important for the en- 
hancement of the ferroelectricity due to orbital hybridiza- 
tion. On the other hand, Moskvin et alM- argued that 
the cluster model in Ref4 is oversimplified to account for 
LiCu 2 2 and LiCuV04 in which the realistic geometry 
configuration of the oxygen atoms was not considered. 
All of these facts motivate the present work to explore 
the effect of bond-distortion on the electric polarization 
in the spiral multiferroic magnets in the framework of 
spin-current model. 

The organization of this paper is as follows. In Section 
HI! we perform quantum chemistry analysis on relevant d- 
orbitals and construct cluster models with two M atoms 
and bridge oxygen atoms. Then we calculate lowest lying 
eigenstates for these cluster models. In Section IIII1 we 
calculate the electric polarization for different clusters. In 
Section llVl we discuss three kinds of chain models formed 
by different clusters and possible ordering of electric po- 
larization. In section [Vj we discuss the relation between 
lattice-driven polarization and magnetic-driven polariza- 
tion and apply our theory to multiferroic copper oxides. 
Section [VTl is devoted to summary. 
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FIG. 1. Relevant electronic states at transition metal atoms. 
Five fold degenerate d-orbitals will split into two fold e g and 
three fold £29 levels. The spin-orbit interaction A splits £29 
levels to a doublet T7 and a quartet T$. In the presence of 
effective exchange U between F7 electrons and local magnetic 
moment, T7 will be further split into spin parallel state \P) 
and antiparallel state |AP). 



and Hund's coupling energy^—. Then Hjj will further 
split the T-? doublet into spin parallel state \P) and anti- 
parallel state \AP), 



IP,-' 



\AP,\ 



^|a)+e*cos^|6) 



cos -z- \a) — e 



i<t>j 



sin -Or I b) 



(4a) 
(4b) 



II. CLUSTER MODELS 

We begin with quantum chemistry consideration. As 
is well known, when the 3rf transition metal atoms are 
placed in an octahedral crystal field, the five-fold degen- 
erate d levels will split into e g and ti g levels. Incorpo- 
rating the on-site spin-orbit interaction A, the e g orbitals 
will not be influenced but the ti g manifold will further 
split into a doublet I\ and a quartet T 8 (see FigfT]). For 
simplicity and following Katsura et al~, we will adopt 
the following T? states as the ground state manifold and 
drop the quartet, 



l a ) = 75 (\dx V ,t) + \d yz ,i) + i \d zx ,i)) , 

\ h ) = 75 i\ d xy,i) - \dyz,t) + i \dzx,t)) , 



(2) 



by assuming the spin-orbit coupling A is the largest en- 
ergy scale of the problem. 

As in Ref4, we further invoke the effective exchange 
interaction between Ty electrons and local magnetic mo- 
ment governed by the following Hamiltonian, 



H v = -Uj^ej ■ S j: 



where e» 



(COS (j)j 



sin( 



:/' 



sin 



vj sin 9j , cos 9j , 



(3) 



is the unit 



vector of the local magnetic moment from the transition 
metal atom M at j-th site, Sj the total spin operator 
of the d-orbital electrons, U the effective exchange in- 
teraction which is of the order of Coulomb interaction 



where \Pj) and \APj) indicate the spin state parallel and 
anti-parallel to the unit vector e} , and the corresponding 
eigenvalues are E\ P .\ = — ¥ and E 
tively. Note that w 
terms of t2 g states, 



tively. Note that we may further write | ± 3 



w v = f, respec- 
P,) and \APj) in 



|P,) =J2a»jm 



<U)\ 

"ftcr/ 



and 



\AP J )=J2 B U)\d$), 



(5a) 



(5b) 



1 1 a 



where /i = xy^yz^zx, a =t;l, A?Z and B^Z are coef- 
ficients obtained from combining Eqs.((4]) and ([5]) which 
depend on angles 6j and <pj (See Appendix [Al for details). 
Hereafter, we assume that A and U are much larger than 
other relevant energy scales so that there is only one rel- 
evant state in each M atoms, say, \P). We also assume 
that there is a spiral magnetic order relating to the tran- 
sition metal atoms. 

Now we are in the position to construct cluster models 
by including bridge oxygen atoms between two M atoms. 
Taking the hole picture, we assume that the oxygen p- 
orbital has energy level E p which is above E\ P \ with a 
difference A as shown in Figl^J i.e. E p = E\p\ + A. The 
hybridization between spin parallel state \P) and oxygen 
p-orbitals results in bonding and anti-bonding states with 
lower and higher energy respectively. In such a M-O- 
M cluster (see FigfJ]), left and right bonding states will 



further hybridize with each other and lead to two lowest 
energy levels Ey and E 2 . When the Fermi energy Ep 
locates between E\ and E 2 , it is called "single hole" in 
Ref4. Whereas, it is called "double hole" when Ep is 
above E 2 but below other levels. 
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FIG. 2. Energy levels demonstration for a M-O-M cluster. 

The hybridization between left and right bonding 
states will also lead to distorted electronic density and 
thereby possible electric polarization in such a M-O-M 
cluster. We may further consider more bridge oxygen 
atoms between two M atoms and derive similar elec- 
tronic density distortion and electric polarization. We 
shall study three-atom M-O-M cluster and four-atom M- 
O2-M cluster in the remained part in this section and 
calculate electric polarization for different kinds of chains 
formed by three-atom or four-atom clusters in next sec- 
tion. 



A. Three-atom M-O-M cluster 

Firstly we consider the three-atom M-O-M cluster 
model shown in Fig|3] Mi and M2 refer to two tran- 
sition metal atoms and the intermediate oxygen atom 
(O) deviates from its centrosymmetric position on Mi- 
M2 bond. To be simple, we assume that the three atoms 
of the cluster are restricted within the xy-plane, the Mi- 
M2 bond is along the £-axis, and the lengths of two M-0 
bonds are equal to each other. Note that Katsura et al£- 
considered a centrosymmetric three-atom model in which 
the oxygen atom locates at the mid point between two 
M atoms. 

Upon the above quantum chemistry analysis, we use 
the following Hamiltonian to describe the on-site energy 
for oxygen p-orbitals and spin parallel states, 

Ha= Yl E \P) c \p d ) c \Pi) +1 ^2( E \P) +A )pl,TP^' ( 6 ) 

.3=1,2 m.ct 

where p^a, (p = X,y,z) are annihilation operators for 
oxygen p-electrons, a is spin index, c\p-\ is the annihila- 
tion operator for spin parallel state at site j. 

The hybridization between p-orbitals and spin parallel 
states is governed by the hopping Hamiltonian, 
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FIG. 3. Three atom M-O-M cluster model. Mi and M 2 de- 
note two transition metal ions with noncollinear magnetic mo- 
ments. O represents a bridge oxygen atom which deviates 
from its centrosymmetric position (denoted by dash circle) 
with a bond-bending angle ip. The arrows on Mi and M2 in- 
dicate the spin direction. The three atoms of the cluster are 
assumed lying in the xy plane. 



where H\ m and 7J 2 m indicate the hopping processes 



between the oxygen atom and Mi and M2 respectively. 
We then write down H\~ m and H^~ m in terms of non- 
vanishing hopping integrals explicitly, 

tfi- m = V(ViP f d (1) +V 2 p^ d (1) 

■ L± t / j \ y *-i J x,a KJj xy,<7 * v ■£t J y,(j^xy,(T 

a 

+ v3pIAI„+v*pI,AI°)i ( g ) 

H 2- m = £ (y lP l a d^ - V,p\A 

a 

-V 3 p{ a dW a + V 4 pi a d^ a ), (9) 



H t = H 



Hi 



H.c. 



(7) 



where d^ a (/i — xy,yz, zx) are annihilation operators 

for d electrons of the transition metal atoms. The su- 

(i) 
perscript j = 1,2 in d]i' a denote the site of magnetic 

transition metal atoms as shown in Fig|3J The integrals 

V n (n = 1,2,3,4) between d and p orbitals are given 

according to Slater-Koster's rules^ as follows, 

Vi = sin (p ( V3 cos 2 ipt pdlJ + (1 - 2 cos 2 ip)t pd A , 

V 2 = cos ip ( V3 sin 2 ipt pdlJ + (1 - 2 sin 2 p)t pd A , 

V 3 = cos (pt pdjr , 

V4 = ahxcptpdw (10) 

Here ip is the angle between M1-M2 bond (i.e., along 
x-axis) and Mi( 2 )-0 bond, as shown in Fig. [3J t pda 
and t pdn are hybridization integrals corresponding to a 
and it bonding between p and d orbitals respectively. In 
contrast to the hopping Hamiltonian derived by Katsura 
et al. in Ref^ where only t pd7T is present due to inver- 
sion symmetry, there are four relevant and independent 
hopping integrals V ni n — 1,2,3,4. Because more ac- 
tive orbitals are involved in our model due to the bond- 
distortion, such &sp x and d yz orbitals. Thus, the Hilbert 
space in the present model is enlarged. The basis of 
the our cluster system contains \Pj), (j = 1,2) and p^ a 
(/i = x, y, z, a —^ , l) and an eigenstate is a linear combi- 
nation of these states. We can easily recover centrosym- 



metric cluster model in Ref4 by setting ip — 0, which The on-site energy in Eq.© is then modified as 

leads to V 1 = V 4 = and V 2 ^V^. 

We proceed to calculate the eigenstates of the system 
by assuming V ^C A and treating H t as a perturbation to 
H a- To the second order perturbation^, the two lowest 
lying states are given by energy levels 



H A = £ E 1p) c\ p . } c IPj) + J2( e \p)+ a )pI ( M^ ( 15 ) 

j — 1.2 i-if-LiC 



E m =2(CT\B\), 



(11) 



with corresponding eigenvectors |1) and |2) given in Ap- 
pendix [B] where the parameters C and B are as follows, 



V? + V 2 2 + Vf + Vi 



C 3A 

B = A ia + A 2 /3, 



(12) 



with 



Ai = - 



v?-v?~vi + vi 



3A 



A 2 = 



ZV3V4 
3A 



(13) 



and 



a = sin — sin h cos — cos — e yvl V2 ' 

2 2 2 2 

= i [sin — sin cos — cos — e Wl w 

M L 2 2 2 2 J 



(14) 



The angles (6i,<fri) and (0 2 ,4> 2 ) specify the direction of 
local magnetic moments in transition metal atom Mi and 
M2 respectively. For a centrosymmetric cluster with ip = 
0, we have A 2 — and C = —A\. 



B. Four-atom M-O2-M cluster 

Secondly, we consider the four-atom M-O2-M cluster 
illustrated in FigJU In this situation, there are two oxy- 
gen atoms symmetrically deviated away from the M1-M2 
bond with bending angle if. Here, we also assume that 
the atoms of the cluster are restricted within the xy plane 
and the M1-M2 bond is along the x-axis for convenience. 



where the superscript i = 1, 2 hxp^ a , (fi = x, y, z) denote 
the site number for oxygen atoms Oi and O2 shown in 
FigJU Note that we assume that p orbitals in the two 
oxygen atoms Oi and O2 has the same energy levels, say, 
E\ P \ + A, due to the symmetric deviation 



E — E 



of two oxygen atoms. 

Accordingly, we shall replace hopping terms H t ~ m and 
H?- m in Eq. © and Eq. © by 



H l-rn = Y ( Kpt(l) d (l) + T/ 2 pt(l) d (l) + Vp Hl) d W 
a 

+ V^Jd%-V iP f)d%), (16) 

and 

H 2-rn = Y (y 1 t(l) d (2) _ ytfMtf*) _ ViP™ <0 



+ v^Jdfl - v lP ^Sil - v 2 p\y d fl a 

(17) 



y ol J z,(T ^zx^a V/ ±l J Z*(T yz, 



2,fT yz,<7 

respectively, where V n (n — 1,2,3,4) are the same hy- 
bridization integrals defined in Eqs. (|TU1) . As shown in 
the previous subsection, due to the oxygen atoms dis- 
placements, there are more active orbitals involved. Rele- 
vant states in such a four-atom cluster include \Pj), \APj) 

and pjlja, where i,j = 1,2, p, = x, y, z, a =t,i- 

In the similar way as the three-atom cluster, the two 
lowest lying states of the four-atom cluster can be ob- 
tained with energy levels, 



E 1(2) =2(CT\A ia \), 



(18) 



and corresponding eigenvectors |1) and |2) given in Ap- 
pendix [C] where the parameters C, A\ and a are those 
given in Eqs. (JI5J), dSJ) and flU). 
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FIG. 4. Four-atom Mi-Oi(2)-M2 cluster. Two bridge oxygen 
atoms Oi and O2 symmetrically deviates from M1-M2 bond 
with a bond-distortion angle <p. The arrows on Mi and M2 
indicate the directions of spins. The four atoms of this cluster 
lie within the xy plane and M1-M2 bond lies along the rr-axis. 



III. ELECTRIC POLARIZATION IN CLUSTERS 

We shall calculate electric polarizations P = (er) for 
three-atom M-O-M clusters or four-atom M-O2-M clus- 
ters using the states |1) and |2) obtained from cluster 
models in previous section. The position of Fermi en- 
ergy Ep is crucial to determine the expectation value P, 
which is given by 

P = n F{El -E^^ + n F{ E 2 -E F )<^, (19) 

where n F (E) = tjtwe is the Fermi function. At low 
temperature, one may replace n F (E) by 6(—E). In the 



case of single-hole, say, E\ < Ep < E 2 , only the lowest 
state |1) contributes to electric polarization. While in 
the double- hole situation, Ep > E 2l both state |1) and 
1 2) are active^. 

According to the expression of 1 1) and 1 2) given in Ap- 
pendix iBl and Appendix [Cl the electric polarization P can 
be written in terms of the overlap dipole matrix elements 
I^ u (a) given as follows, 



I«(a) = e(d^\a\p„(?+a)) 



(20) 



where a is the displacement from a transition metal atom 
M to one of its neighboring oxygen atoms O, a — x, y, z is 
one of the three components of r , /1 = xy, yz, zx denotes 
three t 2g orbitals and v = x,y,z denotes three oxygen p- 
orbitals. We calculate I^ V (S) in Appendix IDl and study 
the configurations of electric polarization in the following 
subsections. 

We find that the electric polarization consists of two 
parts. The first part can be written in terms of ej x e*/+i, 
which is called "spin-current polarization" . The second 
part depends on the bond bending angle and distortion 
configuration, which is called "lattice-mediated polariza- 
tion" . 



A. Three-atom M-O-M cluster 

For a three-atom M-O-M cluster, both spin-current po- 
larization and lattice-mediated polarization contribute to 
the electric polarization. However, they may play differ- 
ent roles in different situations. We shall discuss singe- 
hole and double-hole situations respectively. 

For single-hole situation when E\ < Ep < E%, only 
the lowest state |1) contributes to electric polarization. 
The induced electric polarization Pj.j+i at each element 
cluster connecting j and j + 1-th transition metal atoms 
can be calculated with the help of Eq. (jBll) . The result- 
ing polarization with three components reads, 



px _ n 

py _1 P 1 A 1 +I%A 



(2j,i+i x (e*j x e J+1 ))\ 



"+ 1 ~ 3 A \B jJ+l \ v " +i y 3 3+un y 

2 If 2 J 2 M 1 |q J -, J+1 | 2 +J 1 M2|/3 J -, J+ i| 2 

3 A 3 A\B jd+1 \ 



1 ilAx 

I I!A 2 



P W'+1 = i A 7u L . I @M+1 X (^ X S 3+i ] 



3 A \B Jt j + i 



+ - 



■ (sin 9j cos 9j+i sin < 



3 A \ B j,j+i 

+ cos9jSm0j + ism(j)j + i), (21) 

where -Bj^+i = AiOjj+i + A20jj+i, and ct/j+i and 



Pjj+i are given as follows, 



sin -r- sin ' 



2 
■I ■ e J 



— sin — — 
2 2 



Oj "j+1 -ii 

cos -f cos - J —— e ( 
2 2 

6j 8j + i _ 
cos — cos — — e 
2 2 



-<Pi + i)\ 

(22) 

by replacing cj>i and 4> 2 in Eq. ([T4[) by 4>j and 4>j+i respec- 
tively. And other parameters are defined as 



n = v 3 b 



yz,z 



J 2 



-VI y 

v a xy,; 



+ Volt, 



v 3 iy 



vjy 



yz,zi 



3 1 xy,x ' 2 xy,y ' 3 *zx,. 

I\ = V 2 Iy Z ^y 



*4-lyz,Zl 



VaI; 



xy.z 



Vil 



yz,x- 



(23) 



Each component in the electric polarization Pjj+i con- 
tains two parts, one part can be written in terms of 
Bj x e}+i which has been predicted in Ref£, the other 
part is related to the bond distortion and depends on 
the bending angle if. We denote the second part subject 
to bond distortion as "lattice-mediated", to distinguish 
it from the "spin-current" polarization in the first part. 
Note that, both lattice-mediated contribution and spin- 
current contribution are of the order of V/A and com- 
parable to each other when if is large enough. To derive 
PJ a + i, we use the condition that the two M atoms lo- 
cate mirror-symmetrically about the yz plane across the 
oxygen atom and the T7 manifold in Eq.® which keeps 
unchanged under the mirror operation x — > — x. 

For double-hole situation when Ep > E 2 , two low- 
est states |1) and |2) are involved. Using Eq. (TK)l) and 
Eqs. (|Bl[) (|B2[) . we can easily calculate the electric polar- 
ization Pj j+i, which results in 



Pf 



p y 



0, 



2 If A x + II A 2 



4IJ 

3 A 



C* 



(^) 



(e ij+ i x (ej x e j+1 ))\ y 
4,I^A 1 \a jJ+1 \ 2 + I^A 2 \p jtj+l \ 



P 



2 II A x 

; ■ ' ' ~ 3 A 2 
2 I\A 2 



3 A 2 

(ejj+i x (ej x e j+ i))\ z 



— 2 ' (sin 9j cos 9j+i sin <pj 



i j+1 smcj) j+lj 



(24) 



To derive P, 
-1 



j,j+i-, we have made the approximation 
(111)" 1 --. i + (C-|S|)/A and (2I2)" 1 « 1 + 
(C+\B\)/A. Similar to the single-hole situation, the 
electric polarization can be divided into two parts, the 
lattice mediated contribution and the spin-current con- 
tribution too. However, it is different from the single-hole 
situation that the spin-current contribution is of the or- 
der of (V/A) while the lattice mediated part is of the 
order of V/A. Because the electric polarization contains 
two parts from states |1) and |2) respectively which tend 
to cancel each other—. So the lattice mediated contribu- 
tion will dominate over spin-current contribution when if 
is large enough. 



B. Four-atom M-O2-M cluster 

For a four-atom M-O2-M cluster, only the spin-current 
polarization contributes to the electric polarization while 
lattice-mediated polarization vanishes. Using the results 
in Appendix [C] one can calculate electric polarization 
straightforward. 

For single-hole situation, one obtains that 



P x — fl 



pv ._ 

«' +1 3 A |4i 



jv a x (^J +1 x & x e*j+i)J 



P 



21% A 



A 



\<Xj,j+l\ 

e jtj+1 x (ej x e j+1 )j 



3,3+1 



3AUi 



*J,3+1| 



(25) 



where (Xjj+i is given in Eq. (|22[) . 

For double-hole situation, one obtains that 



^+i=0, 



pv _ 



P 



i,i+i 






(26) 



One sees that lattice-mediated contribution to the elec- 
tric polarization is canceled due to the reflection sym- 
metry between two oxygen atoms in such a M-O^j-M 
cluster. 



IV. CHAIN MODELS AND ELECTRIC ORDERS 

We go further to investigate different chains formed by 
three-atom M-O-M clusters or four-atom M-O2-M clus- 
ters and study possible electric orders in such chains. 
Since most multifcrroics are insulators, relevant electrons 
are localized within clusters, we may neglect electron 
hopping between clusters and use relevant local states 
|1) and 1 2) obtained from cluster models in previous sec- 
tion to calculate the electric polarization P within each 
cluster. 



A. Uniform bending M-O-M chain 

Firstly, we consider the uniform bond-distortion chain 
model shown in FigjS] In this case, all the bridge oxygen 
atoms deviate from the spin chain formed by transition 
metal atoms in the same direction and with the same 
bond-bending angle (p. This type of bond-distortion may 
be ascribed to the DM interaction relating to spiral spin 
ordering, where the spin-helicity vector has the same sign 
for all pairs of neighbouring spins in case of transverse- 
spiral ordering, the DM coupling pushes negative oxygen 
ions in the same direction 5 . 
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FIG. 5. (Color online) A sketch of the uniform bending M-O- 
M chain. All the bridge oxygen atoms deviate from the tran- 
sition metal atoms chain in the same direction and lie within 
the xy-p\a,ne. And the spin-spiral plane is :ry-plane. The spin- 
current and lattice-mediated polarization are denoted by red 
and green arrows, respectively. 



It is easy to see that ferroelectric order will form in 
this kind of chains. To illustrate it, we assume all the 
magnetic moment from transition metal ions lie in xy- 
plane and are spiral ordered in the same plane as shown 



in FigEl In this case, Pf,- + i = 
does not depend on site index j. 
For single-hole, 
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while for double-hole, 
2 If Ax - 
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with A(f> — <j)j ~ 



>j+i- 



B. Staggered bending M-O-M chain 

There also exist some materials with staggered bend- 
ing M-O-M chains as shown in Fig|6] For instance, the 
Mn-O-Mn-O-Mn bonding along (110) direction forms a 
zigzag chain due to alternative rotation and tilt of the 
MnOe octahedra in some perovskite rare earth man- 
ganese oxides RMnO^i. 

In this case, spin-current polarization and lattice- 
mediated polarization may form different orders. To il- 
lustrate it, we still assume all the magnetic moment from 
transition metal ions lie in xy-plane and are spiral or- 
dered as shown in Fig|6l In this case, Pfj+i = Pjj+i — 
0. Pjj+i contains two parts, spin-current polarization 
Pj C j+i an d lattice- mediated polarization P,7+i' 



pV _ psc,y 



P 



lm,y 



3,3 + 1- 



(29) 
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FIG. 6. (Color online) A sketch of the staggered M-O-M 
chain. The neighbor oxygen atoms of a given magnetic metal 
atom displace toward opposite directions. The spin-current 
polarization is uniform in each M-O-M cluster (red arrows), 
while lattice-mediated polarization is staggered (green ar- 
rows). 



For single-hole, 
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while for double-hole, 

2i 1 , A 1 +i|A 2 . 
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So that spin-current polarization is uniform and forms 
ferroelectric order while lattice-mediated polarization is 
staggered and forms antiferroelectric order. The total 
electric polarization is "ferrielectric ordered". 



C. Symmetric M-0 2 -M chain 

Finally, we consider a symmetric M-O2-M chain shown 
in FigtZl which is made of four-atom M-O2-M clus- 
ters. Such symmetric M-O2-M chains exist in two 
prominent examples of multiferroic compound, namely, 
LiCu202^^ and LiCuVO^, which crystallize as one- 
dimensional chains of the edge-sharing CUO4 plaquettes. 

In this kind of chains, lattice-mediated polarization 
vanishes and spin-current polarization may form ferro- 
electric order. Assuming all the magnetic moment from 



FIG. 7. (Color online) A schematic edge-shared MO2 chain. 
The bridge oxygen atoms are symmetrically located beside 
the two sides of the magnetic metal atoms chain. In this case, 
the lattice-mediated contribution of the electric polarization 
is exactly canceled. The nonvanish part of the polarization 
arise from the spin current, indicated by red arrows. 



transition metal ions lie in xy-plane and are spiral or- 
dered as shown in FigfTl we obtain Pfj+i = Pjj+i = 
and 



_ 2 If A, sin(A0) 
«' +1 3 A |Ai 



cos • 



in single-hole situation and 



■\w-t-i ~ 
in double-hole situation. 



8 IJM 
3 A 2 



sin(A0), 



(32) 



(33) 



V. DISCUSSIONS 

In this section, we shall discuss the relation between 
magnetism-driven electric polarization and lattice-driven 
electric polarization, and then apply our theory to mul- 
tiferroic copper oxides. 

Lattice-driven electric polarization: As shown in FigJSJ 
the negatively charged oxygens coherently displace away 
from the magnetic chain formed by M ions. This shift 
will generate a net lattice- driven electric polarization 
pid ducting from the oxygen ions to the M ions chain. 
The polarization P ld can be evaluated as long as the 
p-d hybridizations t p da and t p dn and the geometry of 
the lattice are specified. We shall compare this lattice- 
driven electric polarization (P ld ) with the two parts of 
magnetism-driven electric polarization studied in this pa- 
per, say, lattice-mediated polarization (P Jm ) and spin- 
current polarization (P sc ). To do this, we restrict our- 
selves to the double-hole situation and choose the follow- 
ing parameters, tpda/tpdir = —2^ A = 2eV, sin(A</>) = 
sin(0.287r)i2ii2iI4 and the lattice constant a = 5k. The 
results are shown in Fig|8] where electric polarizations 
are plotted as functions of bond-bending angle ip. For 
a typical ip, P ld ~ 10 2 ^C/m 2 , while P lm ~ 10 2 /iC/m 2 
and P sc ~ 10/iC/m 2 . The lattice- mediated polarization 
keeps in parallel with the lattice-driven polarization for 
any bond-bending angle p while the spin current polar- 
ization may varies with tp and be in parallel or anti- 
parallel to the lattice-driven polarization. Moreover, the 
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spin current polarization depends on spin configuration 
and can be neither parallel nor anti-parallel to lattice- 
driven polarization. Upon these observations, we con- 
clude that lattice-driven polarization is cooperative with 
lattice-mediated polarization while spin current polariza- 
tion may be either competitive or cooperative, depending 
on the spin configuration and the bond-bending angle. 

Note that lattice-driven polarization P ld in FigJS^a) is 
overestimated since we chose the maximum electric po- 
larization in all the possible valences for magnetic transi- 
tion metal and oxygen ions. It is also worth noting that 
altough P ld and P lm dominate over P 3C at finite ip in 
uniform bending chains as shown in Fig[8l P sc will be 
dominating at small ip or in symmetric M-O2-M chains. 
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FIG. 8. (Color online) The ^-dependence of electric polar- 
izations in uniform bending M-O-M chains (see Fig[5}. (a) 
Lattice-driven electric polarization, (b) Lattice-mediated po- 
larization, (c) Spin-current polarization. 

Multiferroic copper oxides: Then we shall apply our 
theory to several copper oxides, including LiCuVO,!, 
LiCu2C>2 and NaCu2C>2, where magnetic copper ions 
Cu 2+ and oxygen ions O 2 " form symmetric M-O2-M 
chain as shown in FigJT] 

There are two prominent compounds which can be 
described symmetric M-O2-M chains, LiCu2C>2 and 
LiCuVC>4. Both of them are characterized by edge- 
sharing CuC>4 plaqucttes, forming spiral spin orders and 
emerging electric polarization at low temperature. The 
electric polarization in these two compounds is rela- 
tively weak comparing with manganites, ~ 10/iC/m 2 , 
and shows strong anisotropy. There is another interest- 
ing material NaCu202, which is isostructural to LiCu202 



while Li is substituted by Na. Although NaCu2C>2 also 
exhibits helical magnetic order at low temperature^—, 
no electric polarization has been observed in it within the 
experimental sensitivity (<0.3/iC/m 2 ) 39 ' 40 . We shall ap- 
ply our theory to these prototypical compounds and ad- 
dress the issue of "missing multiferroicity" in NaCu202 
in the rest of this section. 

As mentioned in Section IIIIB1 the lattice-mediated 
contribution is exactly cancelled in a symmetric M-O2- 
M cluster, and the electric polarization comes from the 
spin current only. In contrast to widely used isotropic 
form of Eq.([T|), the spin current form of electric polar- 
ization in the symmetric M-O2-M cluster is anisotropic 
in general, which is in accord with the experimental 
finding o 24 ' 25 ! 27 ' 28 and density functional calculations 41 . 
The magnitude of the electric polarization depends on 
the bond-bending angle tp and vanishes at (p = vt/4. Note 
that the spin current form will recover isotropic in the 
absence of bond distortion, say, tp = 0. 
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FIG. 9. The (^-dependence of y- and z- components 

of electric polarization in the presence of xy- and 2x-plane 
spin spiral order respectively in symmetric M-O2-M chains 
(see Fig[7]). A typical set of parameters are choosen for 
Cu 2+ oxides^, t pdn = 0.75eV, t pdtT = -1.85eV, A = 3eV. 
The Clementi-Raimondi effective charges Zq = 4.45 and 
Zq u = 13.2 are used for O 2- and Cu 2+ ions respectively^ 
and Rcu-o = 1.445A (i.e., LiCuVO^). 

In FiglHl we plot y- and z- components of electric polar- 
ization as functions of bond-bending angle tp in the pres- 
ence of xy- and zz-plane spin spiral order respectively. 
One sees that the polarization is sensitive to the bond 
bending angle <p and is anisotropic. In realistic materials, 
the bond angles of Cu-O-Cu in LiCuV04, LiCu2 02, and 
NaCu 2 2 are 95°^, 94°^, and 92.9°^^ respectively, 
corresponding to bond bending angles ip = 42.5°, 43°, 
and 43.6°, which are all close to ip = tt/4. As shown in 
the Fig. [9] the electric polarization decreases and ulti- 
mately vanishes when ip approaches vt/4. Indeed, exper- 
iments clearly indicated that the electric polarization in 
LiCuVC>4 is stronger than that of LiCu2 02- Moreover, 
the bond-bending angle <p in NaCu202 is closer to tt/4 
than that of LiCuV04 and LiCu202, resulting in weaker 
electric polarization which is difficult to be observed. 

For comparison, we also list calculated electric po- 



TABLE I. The a- and c- components of electric polarizations 
in LiCuVCU, LiCu2 02 and NaCu2C>2 in the presence of ab- 
and fee-plane spin spiral orders. The values in parentheses 
correspond to experimental data. 



Materials 


<P 


Polarization (/iC/m 2 ) 
P a (ab) P c {bc) 


LiCuVO^^ 
LiCu 2 02^^ 
NaCuaO^ 


42.5° 

43° 

43.6° 


5.8 (43) 4.3 (10) 
4.4 (8) 3.8 (4) 
2.1 (?2) 0.8 (0) 



[^Currently no experimental data is available. 



larizations for these three compounds in Table U using 
their concrete parameters. For example, the pitch an- 
gles for LiCuV0 4 , LiCu 2 2 and NaCu 2 2 are 84.2°^L2£, 
62.6°—, and 81.7°— respectively. Corresponding Rrj u -o 
(Rcu-o is the distance of Cu-0 bond along the Cu-Cu 
chain) are 1.445AM, 1.43A2£, and 1.465A 4 -2 respectively. 
Other parameters are the same as those in FigEI The 
values in parentheses in Table U are experimental data. 
The decreasing electric polarization with ip approaching 
7r/4 is in agreement with experimental observation. For 
NaCu 2 2 compound in the presence of fee-plane spin spi- 
ral orde r 36 ' 38 , we find that P c is extremely weak and hard 
to detect^. To our knowledge, there is no experimental 
report on electric polarization in the presence of afe-plane 
spin spiral order, which may be tunable under external 
field. Our theory predict a finite electric polarization 
(~ 2/j,C/m 2 ) oriented along the a-axis in the presence of 
a6-plane spin spiral order. This prediction can be tested 
in future experiment for NaCu 2 2 . Thus, the effect of 
bond-distortion on the magnetic-driven polarization pro- 
vides a clue to clarify the anisotropy of the electric polar- 
ization observed in LiCuV04 and LiCu 2 2 and to eluci- 
date the puzzle of non-multiferroicity in NaCu 2 2 . 



distortion is present. For four-atom M-0 2 -M clusters, 
lattice-mediated polarization vanishes due to symmetry 
and there exists only spin-current polarization which is 
of the order of (V/A) in single-hole situation and of the 
order of (V/A) 3 in double- hole situation. 

Then we study three kinds of chain models made of 
different clusters where magnetic moments are spiral or- 
dered. Firstly, we study uniform bond-distortion chain 
models where all the bridge oxygen atoms deviate from 
the spin chain formed by transition metal atoms in the 
same direction and with the same bond-bending angle. 
Ferroelectric order may exist in such chains. Secondly, 
we study staggered bending M-O-M chains in which the 
bridge oxygen atoms displace on opposite directions be- 
tween two adjacent clusters. Ferrielectric order may ap- 
pear instead of ferroelectric order in this kind of chains. 
Thirdly, we study symmetric M-0 2 -M chains which are 
made of four-atom M-0 2 -M clusters and the two oxygen 
atoms deviate symmetrically from the M-M bond within 
a M-0 2 -M cluster. The electric polarization may be fer- 
roelectric ordered in such symmetric chains. 

We also discuss the relation between magnetic-driven 
electric polarization and lattice-driven electric polariza- 
tion in the presence of bond distortion in details. 

Finally, we apply our theory to multiferroics copper ox- 
ides and find the results agree with experimental observa- 
tions. The issue of "missing multiferroicity" in NaCu 2 2 
compound is clarified and more experimental predictions 
are made. 

Summarizing, we show that the M-O-M bond distor- 
tion may affect electric polarization in multiferroics sig- 
nificantly and the possibilities for ferroelectric and ferri- 
electric ordering are discussed. The application to mul- 
tiferroics copper oxides is successful. 
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Appendix A: Coefficients A% and fife in Eqs.([5]) 



In this appendix we shall derive coefficients A 1 , -, and 
B 1 ,^ in Eqs. (l5al) and ([5b]). For the parallel state \Pj) and 
corresponding coefficients A Z ,' T ) , we combine Eqs. (f4"al) and 
(l5al) to have 



sin -j\a) + e i **cos^\b) = '£ A U) 



<&)) 



(Al) 



Using the relation Eq.Q, we obtain A\ a -\ as follows 
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Similarly, we obtain coefficients B 1 ,?-. for the anti- 
parallel state \APj) as follows, 
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Appendix B: Lowest lying eigenstates for three-atom 
M-O-M cluster 



We calculate the eigenvectors for two lowest lying 
eigenstates |1) and |2) in a three-atom M-O-M cluster 
in this appendix. The Hilbert space contains states \Pj), 
(j = 1,2) and p^ a (/j, = X,y,z, cr =t,|). Assuming 
^ « A and treating H t as a perturbation to Ha, we 
obtain the two eigenvectors |1) and |2) up to the second 
order perturbation, 



B 



V2\B 
1 

with Ei =2(C- |B|), and 
B 



cr 



M^P |p, 



V2|B| 

1 



cr 

i P2 > - 1 E ( A ( 2 T ( Fi I**.*) - y2 1***)) + (-^7 + Msf) b z , 



(Bl) 



(B2) 



with B2 = 2(C + |-B|), where parameters C and B are Appendix C: Lowest lying eigenstates for four-atom 

M-O2-M cluster 



In this appendix, we calculate the eigenvectors for two 
lowest lying eigenstates |1) and |2) in a four-atom M-O2- 
M cluster. In this situation, the Hilbert space contains 



given in Eqs. ([12]) , ([13]) and ([14]) in the main text, coeffi- 
cients A 1 / 7 , are given in Eqs. (|A2p . It is easy to verify that 
the states |1) and |2) are orthogonal to each other and 
their normalization factors are (1|1) = 1 — (C — \B\) /A 
and (2 1 2) = 1 - (C + \B\) /A respectively. 



.,(*) 



states \Pj) (j = 1, 2) and pfc a (i = 1, 2, /i = x, y, z, a =\ eigenvectors |1) and |2) as follows, 
, 4-) ■ Up to the second order of t , we obtain the two 
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with E x = 2(C - |Aia|), and 
j4iq 
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with E% = 2(C + |Aia|), where parameters C, A\ and for nonvanishing I" (a), 
a are given in Eqs. (fT"2|) . (TT3|) and (|14[) in the main 
text, coefficients A l , a \ are given in Eqs. (|A2p . The nor- 
malization factors are (1|1) = 1 — 2 (C — |-4.ia|) /A and 
(2|2) = l-2(C + |i4ia|)/A. 



Appendix D: Overlap dipole matrix elements I^ v (a) 



From the definition of the overlap dipole matrix ele- 
ments I" „ (a) in Eq. (1201) , and using the mirror symmetry 
about the xy-plane (z — ¥ —z), we find that the following 
matrix elements vanish, 
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xy,z 

jy 
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yz,x 

jy 

yz,x 

p 



r 

yz,y 
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zx 
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zx,y 
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zx,y 
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0, 
0, 



xy,x xy,y yz,z 

Then we study symmetry relations for nonvanishing 
I" „ (a) with two vectors a and b, where a and b lie within 
the a;y-plane and are linked by symmetry P x (x — > — x) 
or P y (y — > —y)- When b = P x a or b = P y a, we have 

T x (a) = I x (b) P (a) = -P (b) 

^xy.yy^ 1 / ^xy^yVji^xy^xV* 1 / ^xy^xV^J^ 

Iy Z ,z(a) = -Iy ZiZ (b),I* x Ja) = I* XiZ (b), 

I V xy,x(*) = I V xy,x(b),P xy Jd) = ~P xy Jb), 
I y y Z ,M = I V y Z Ab),P zx Ja) = -P zx Jb), 
KyA®) = - I xy,z( b )Jyz,x(^) = -^,x( & )) 
^«,»(3) = Iy Z ,y(b),Ilx,x( d ) = ^zxjb), 

I z zx, v (a) = -PzxJb), 



It turns out that relevant /" (a)'s in the main text 
include P P P P P P P We de- 

m^iLiia^ ±y Z y, -t^y^, J -y Z ^ x t A X y lX -j ^yZjZi ^zx,z^ A xy,y vv ^ ^^ 

mostrate the numerical evaluation for these integrals in 
Fig nUl and FigfTTI respectively. Rm-o is the M-0 bond 
distance along the M-M chain. In this calculation, we 
have taken the hydrogenlike radial-wave functions and 
the Clementi-Raimondi effective charges Zq H = 4.45 and 
Z$ = 10.53 for 2 ~ and Mn 3+ ions respectively^. It 
is shown that, for typical M-0 separation (i.e., Rm-o ~ 
4ao, here a is Bohr radius), I XVyX is one order of magni- 
tude larger than those of other integrals. The separation 
between the M and O ions Rm-o n ot only yields violent 
variance on the magnitude of the overlap dipole integrals 
but also may change their signs. 
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FIG. 10. Two-site dipole integrals I^ v against M-O separa- 
tion Rm-o with various bond bending angle. The inset in each 
figure shows the zoom in view of the integrals near typical M- 
O separation (i.e., 4ao, io is the Bohr radius). 
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FIG. 11. Two-site dipole integrals I* v against M-O separa- 
tion Rm-o with various bond bending angle. The inset in each 
figure shows the zoom in view of the integrals near typical M- 
O separation (i.e., 4oo). 
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